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Abstract 

We determine conditions for the quantisation of graphs using the Dirac operator 
for both two and four component spinors. According to the Bohigas-Giannoni-Schmit 
conjecture for such systems with time-reversal symmetry the energy level statistics 
are expected, in the semiclassical limit, to correspond to those of random matrices 
from the Gaussian symplectic ensemble. This is confirmed by numerical investigation. 
The scattering matrix used to formulate the quantisation condition is found to be 
independent of the type of spinor. We derive an exact trace formula for the spectrum 
and use this to investigate the form factor in the diagonal approximation. 
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1 Introduction 



Quantum graphs have proved an important model in the semiclassical study of systems 
whose classical analogues are chaotic [T(| |T7|. Much work in this field concerns the con- 
jecture of Bohigas, Giannoni and Schmit || which connects statistics of the energy level 
spectrum to those of random matrices, the ensemble of random matrices depending on 
the symmetries of the system. For systems with no time-reversal symmetry statistics 
of the Gaussian unitary ensemble (GUE) are expected whereas those of the Gaussian 
orthogonal (GOE) or symplectic (GSE) ensemble should apply for systems with time- 
reversal symmetry, depending on whether the spin is integer or half integer, respectively. 
Most investigations of this conjecture have concentrated on systems with spin zero and 
there are few examples of systems whose energy level statistics follow the GSE prediction 



TT|, [L9L pTl O]. In line with this the usual graph quantisation applies the Schrodinger 



operator to a metric graph, finding energy level statistics of the GOE or GUE |16, 17 



It is our aim to develop a quantisation of graphs including half integer spin and we 
expect to observe spectral statistics following the GSE when time-reversal symmetry is 
present. We choose a Dirac operator that we realise as a self-adjoint operator on an 
appropriate Hilbert space. The Dirac operator on a graph was considered previously by 
Bulla and Trenkler [1(J as an alternative model of a simple scattering system, however, 
without addressing the problem of time-reversal invariance. Instead we take closed graphs 
and find boundary conditions that ensure a self-adjoint realisation of the Dirac operator 
such that time-reversal symmetry is preserved. By looking at systems of this type we 
can compare the spectrum with the results for the Schrodinger operator on graphs and 
general semiclassical results for systems with spin. One well known advantage of the 
quantum graph is that the trace formula is exact rather than semiclassical. This provides 
the opportunity to distinguish features present only in the semiclassical limit from those 
inherent in systems with spin. 

One unusual quality of the Dirac operator in one dimension is the possibility of two 
component spinors rather than the usual four component spinors required in three di- 
mensions. Physically this appears odd. If the Dirac operator on the graph describes an 
idealisation of a physical system of wires in the limit that the width of the wire tends to 
zero we are lead to four component spinors of the type required in three dimensions. If 
however we set up the mathematical problem in which the graph is a topological entity 
then it is natural to choose two component spinors. In the context of quantising graphs the 
apparent contradiction will in fact disappear. The spectrum is independent of the choice 
of spinors so the physics of the system cannot distinguish the language used to describe it. 

Our paper is organised as follows: In section ||] we introduce the necessary terminology 
of graphs. Section ^ sets out the two approaches to the Dirac operator in one dimension, 
restricting from the Dirac equation in three dimensions which leads to four component 
spinors or considering an irreducible representation of the Dirac algebra which requires 
only spinors with two components. In section [| we find self-adjoint and time-reversal 
symmetric realisations of the Dirac operator for two component spinors and in section 
[5] we do the same in the four component case. Numerical examples confirm that even 



2 



for simple graphs the energy level statistics follow the GSE prediction. We find that the 
scattering matrix in both formalisms has the same properties, the spectrum is therefore 
independent of the approach taken. In section ^ we first derive the exact trace formula for 
the density of states and then use this to calculate the form factor. Making the diagonal 
approximation we see that the form factor agrees with the GSE result for low r. Finally 
we suggest how the recipe for quantising a graph with the Dirac operator is consistent with 
the model of a graph as a three dimensional system of wires, section |7|. 



2 Graphs 

We begin with a few general definitions necessary when considering graphs. A graph 
consists of V vertices connected by B bonds. The valency v i of a vertex i is the number 
of bonds meeting at the vertex. The topology of a graph can be described using its 
connectivity matrix C. This is a V x V matrix with entries 

q ( 1 if vertices i and j are connected, (2 1) 

13 \ otherwise, 

where we have assumed that the graph has at most one bond connecting any pair of 
vertices. We also suppose the graphs under consideration to be connected, any vertex can 
be reached from any other by passing down bonds. A graph is directed if the bonds are 
assigned directions: A bond running from vertex % to vertex j will be denoted as b = (ij). 

On a graph one usually considers probabilistic rather than deterministic classical dy- 
namics. This can be described by a stochastic matrix M of transition probabilities between 
the bonds. The matrix M = (M&c) defines a Markov chain on the graph propagating vec- 
tors p = (pi, . . . , P2b) T of probabilities assigned to the bonds, p i— > Mp. Note that we allow 
particles to move in either direction on a bond and so there are 2B possible states of the 
particle on the graph. (This convention for counting bonds will prove necessary later as 
plane-wave solutions for the Dirac operator cannot be made to correspond to the directions 
on a directed graph as in the Schrodinger case.) The non-negative matrix M is called irre- 
ducible if for every pair (b, c) of states there exists a power q such that (M q ) bc > 0. If the 
power q can be chosen independent of (b, c), such that all elements of the matrix M q are 
strictly positive, M is called primitive. A Markov chain with an irreducible Markov matrix 
M is ergodic, i.e. any initial distribution p converges in average to the stable distribution 

P° = 2b(1> ' ' ' ' 1) T ' 

N 



& = (2-2) 

n=l 

A Markov matrix corresponding to an ergodic chain has a non-degenerate eigenvalue one. 
For primitive stochastic matrices all other eigenvalues are smaller in magnitude such that 
the associated Markov process is mixing. Then an initial distribution converges to the 
stable distribution, 

lim M n p = p , (2.3) 
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and, moreover, 



lim (M n ) bc = -L . (2.4) 

n— >oo ZlD 



It is this last condition that will be used in section |6.2| to evaluate the form factor in the 
diagonal approximation. Details on non-negative matrices and Markov chains can be found 
in 
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If the quantum system on the graph is defined by a matrix T = (T{, c ) of complex 
transition amplitudes between the bonds (we will see that T is unitary) then 

M (ij)(ki) = \ T (ij)(ki)\ 2 ■ (2.5) 

M(ij)(ki) is the quantum mechanical probability that a particle traveling from k to i scatters 
at i to travel towards the vertex j. The unitarity of T implies that M is doubly stochastic 
and therefore defines a Markov chain. 

While only topological information is necessary to investigate classical dynamics on a 
graph quantum mechanics requires us to assign lengths, to bonds, (ij). This defines 
a metric graph. It is natural to consider a bond with length as directed, one vertex lying 
at zero on the bond and the other at L, consequently our metric graphs are also directed 
graphs. In order to avoid degeneracies in the lengths of periodic orbits on the graph we 
assume the lengths assigned to the bonds are incommensurate, not related by a rational 
number. 



3 The Dirac equation in one dimension 

In one spatial dimension the Dirac equation is 

d ( d \ 

ih—^(x,t) = -ihca — + mc 2 (3 V(x,t) , (3.1) 
at \ ox J 

where a and (3 satisfy the relations a 2 = (3 2 = I and a(3 + (3a = that define the Dirac 
algebra. There are two possible interpretations of this equation: Either one views it as a 
restriction of the Dirac equation in three dimensions, or one considers it as a problem in one 
dimension from the outset. In the first case this equation originates from an implementation 
of the Poincare space-time symmetries in relativistic quantum mechanics such that the 
notions of spin and of anti-particles can be carried over. If then, for instance, one restricts 
to the y axis the Dirac matrices a = a y and (3 are hermitian 4x4 matrices that form a 
reducible representation of the Dirac algebra in one dimension. The second case is void of 
the physical interpretations deriving from 3 + 1-dimensional space-time symmetries. One 
merely considers the Dirac equation (|3.1|) with a faithful irreducible representation of the 
Dirac algebra and thus chooses a and f3 to be hermitian 2x2 matrices. For an extensive 



discussion of the Dirac equation see [22 



The two cases are naturally connected as a unitary transformation of a y and (3 brings 
them into block diagonal form preserving the algebraic relations. For example the standard 
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(Dirac) representation of the matrices in three dimensions, restricted to the y axis, is 
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(3.2) 



Now consider the unitary transformation 
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(3.3) 



Applying U to a and (3 generates a new representation, 



UaU 



-i 



( 


— i 
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-1 
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- 1 / 



(3.4) 



which is obviously reducible, 
in one dimension irreducibly. 



Either of the 2x2 blocks represents the Dirac algebra 
It is, however, important to note that passing from the 
four dimensional to the two dimensional representation the usual interpretations of spin, 
particles and anti-particles, and time- reversal are lost. 

For both two- and four- dimensional representations of the Dirac algebra the Dirac 
operator reads 



V 



•J- U 2 n 

-lhca - — h mc p 

ax 



(3.5) 



On the real line T> is defined as an essentially self-adjoint operator on the Hilbert space 
L 2 (M) <g> C n with domain C^(R) ® C n , where n = 2 or n = 4. 

In order to obtain self-adjoint realisations of the Dirac operator on a compact interval 
/ = [0, L] appropriate boundary conditions have to be specified. These can be classified 
by extensions of the closed symmetric operator T>° given as T> defined on the domain 
W^ x {I)®C n C L 2 (I)®C n . Here W^I) denotes the Sobolev space of functions ip on I that 
are, along with their (generalised) first derivatives, in L 2 (I) and satisfy y(0) = = <p{L). 
The closed symmetric extensions of T>° arise from restrictions of D defined on the domain 
Wi,i{I) <8> C n of n-component Sobolev spinors with no specified boundary conditions to 
subspaces on which the skew-Hermitian quadratic form 



(3.6) 



vanishes. The maximal isotropic subspaces of W2,i(-0 ® C n with respect to Q, i.e. the 
maximal subspaces on which Q vanishes, then yield domains on which D is self-adjoint. 
This is the approach adopted by Kostrykin and Schrader for the Schrodinger operator |T5[ 
on graphs. Self-adjoint realisations of the Dirac operator on an interval are also classified 
using a different technique by Alonso and De Vincenzo [|TJ. 
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4 The Dirac operator for two component spinors 



From the two approaches to the problem we begin with the minimum necessary and con- 
sider two component spinors first. According to the above we take 

-(!?). />-(;_°0- («> 

(The choice of a and (3 is designed to simplify the calculations.) 

Eigenspinors i/j = (^) of V are solutions of the time independent Dirac equation 

Vip{x) = Eip(x) . (4.2) 

For positive energies E they are hence plane waves of the form 



M*) = m ^ h{k) j + m ^ _ h{k) J e— (4.3) 
with k > and 

7 (Jfe) := E ~ mc E = J (tick) 2 + m 2 c 4 . (4.4) 

The choice of the coefficients fi(k) and (x{k) depends on the boundary conditions imposed 
by the self-adjoint realisation of T>. We see immediately that the plane-wave spinor ( |4.3| ) is 
not invariant changing x to —x. As the Dirac operator is first order the direction assigned 
to bonds on the graph becomes significant. 

4.1 Self- adjoint realisations on graphs 

Self-adjoint realisations of the Dirac operator on graphs were considered by Bulla and 
Trenkler ||10|| . However, in order to develop a simple form related to the vertex transition 
matrix, we continue with the approach of Kostrykin and Schrader for the Schrodinger 
operator. 

The basic idea is to view a graph as a collection of intervals that are glued together 
according to the connectivity matrix C. The Hilbert space on which the Dirac operator 
acts is therefore the direct sum of the Hilbert spaces for each bond, 

B 

H:=Q)L 2 ([0,L b ])®C 2 , (4.5) 

6=1 

where b runs over all bonds. The spinors ip G Ji therefore consist of the B components 
. . . ,ip B ) attached to the bonds, each of which is a two-spinor. The scalar product of 
spinors on the graph is then 

B 

(0,^> = ^(^V (4.6) 

6=1 
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where {(fi b ,ip b )b is the scalar product in L 2 ([0,Lb}) <g> C 2 . Self-adjoint realisations of V are 
constructed in analogy to the case of a single interval by suitably extending the closed 
symmetric operator T>° defined on the domain 



B 



0Wji([O,I*])®C a . (4.7) 



6=1 

The domains of its closed symmetric extensions are the isotropic subspaces of 

B 



0^ 2jl ([O,L b ])®C 2 (4.8) 
6=1 

with respect to the skew-Hermitian form Q as given in (|3.6|) . Integrating by parts yields 

B 

0(&VO =hcJ2(ri(U)t(0)-rt(L b )t(L b ) -^(0)^(0) + $(L h )&(L b j) , (4.9) 



showing that Q only depends on the boundary values of the spinor components. Self- 
adjoint extensions again correspond to maximally isotropic subspaces and thus can be 
characterised by boundary conditions. 

Following the construction of Kostrykin and Schrader [IB*] for Schrodinger operators 
on graphs, we introduce a map from the infinite-dimensional Hilbert space Ti, to the 4B- 
dimensional space of boundary values, ip \— > ij) — f ^ ) ' w ^ n 



T 

T 



1>! := W(0), • • • , (0), Vi(^i). • • • > (^b)) 
V> 2 := M(0), • • • , -^(0), ^(£1), • • • , V>f(£fl)) 
By including ■02(0) with a minus sign we can write Q as 



(4.10) 



"<**>=<*! *d(-L l2 o B )U) • (4 - u) 

The maximal isotropic subspaces in (|4.8|) with respect to Q are now equivalent to maximal 
subspaces of vectors -0 G C 4B on which the complex symplectic form fl4.11| ) vanishes. 
Defining a linear subspace of C 4S as the vectors if) satisfying 

A<0! + BV> 2 = , (4.12) 

with complex 2B x 2B matrices A and B, one first observes that in order for this subspace 
to have the desired dimension 2B the 2B x 4B matrix (A, B) must have maximal rank. 
Kostrykin and Schrader show that such a subspace is maximally isotropic if and only if 
(A, B) has maximal rank and AB^ is hermitian. 
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Equation ( |4.12| ) defines boundary conditions for the Dirac operator on a collection of 
B bonds. These boundary conditions yield a self-adjoint realisation of T> when 

rank(A,B) = 2B and AB f = BA t . (4.13) 

So far the connectivity of the graph has played no role. In order to implement this we 
consider only boundary conditions, specified through the matrices A and B, that connect 
bonds according to the connectivity matrix C. The matrices A and B then have a block 
structure, each block defining boundary conditions at a single vertex. 




Figure 1: An isolated vertex on which we wish to define the transition matrix. 



At a given vertex i we label the Vi bonds such that the first Wi bonds are outgoing and 



the remainder incoming, see figure |l| The vector 
at i then consists of the components 



of boundary values 



Wi+l): 



(-^(0), . . . , -V^(0), i>T + \L m+1 ), . . . , ^(Aj) . 



(4.14) 



If we denote the v j x Vi blocks of the matrices A and B that define the boundary conditions 
at i by Aw and BW, we have the condition 

A^V? + B^V? = with A w B wt = B (i) A wt (4.15) 

which guarantees current conservation at the vertex %. 

Eigenspinors ip of the Dirac operator T> with given boundary conditions are composed 
of plane waves, equation (|4.3|) , on each of the bonds. Assigning incoming and outgoing 
plane- wave solutions to each of the bonds we can write down vectors of coefficients ~jx , *fl 
for the outgoing and incoming waves at the vertex, respectively: 



/ 



-1 g-ifcLt, 



\ p, Vi e 



( 



\ 



/i 



\ n Vi e 



(4.16) 
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Here we omit the index i labeling the vertex. Writing tp\ and ij) 2 in terms of /x and fi 
and substituting into the definition of the boundary conditions we obtain 

Jt = -(A (i) + i 7 (A;)B (i) )- 1 (A (i) - i 7 (A;)B (i) ) ^ . (4.17) 

Notice that under the conditions given the matrices (A^ ± i'j(k)M^) are always invert- 
ible ||15|| . This defines the vertex transition matrix TW which connects the coefficients of 
incoming and outgoing plane waves at the vertex i, 

T (i) := -(A w + i 7 (A;)B {i) )- 1 (A w - i 7 (A;)B (i) ) . (4.18) 
The condition that A^EW i s hermitian implies that the matrix TW is unitary, 

( T «)-i = ( T «)t _ ( 4 _ 19 ) 

Self-adjoint realisations of the Dirac operator on the graph are defined by matrices A, B 
with AB^ hermitian. These prescribe the boundary conditions ( |4.12|) and determine unitary 
vertex transition matrices (|4.18[ ) for plane waves on the graph. 



4.2 Time-reversal symmetry 

Following the conjecture of Bohigas, Giannoni and Schmit, we expect to find energy level 
statistics corresponding to those of the Gaussian symplectic ensemble (GSE) for quantum 
systems with half-integer spin s and time-reversal symmetry. This is due to the fact that 
the time-reversal operator is anti- unitary and squares to (— l) 2s , imprinting a symplectic 
symmetry on the Hamiltonian. In the case of a two component Dirac equation in one 
dimension the usual interpretations of spin and time-reversal inherited from three dimen- 
sions are lost. Nevertheless, one can introduce an anti- unitary operator squaring to — 1 
that is formally identical to the time-reversal operator for two component (Pauli-) spinors 
in three dimensions, 

T ■■=(-! J)*' ( 42 °) 

where /C is complex conjugation. This can be interpreted as a physical time-reversal in 
the sense that it changes the sign of time and momentum. The matrix part of T is only 
required in order that T 2 = —I as there is no spin operator for two component spinors. 
In the subsequent work we will therefore refer to ( f4.20|) as a time-reversal operator. A full 



discussion of time-reversal symmetry with spin is found in the second chapter of 

We now determine conditions under which the vertex transition matrix is time-reversal 
symmetric. This will raise questions about when general time-reversal symmetric boundary 
conditions exist for a Dirac operator on a graph. For the system to be time-reversal 
symmetric T must commute with the Hamiltonian T>. This forces the mass term in D to 
vanish. From now on we consider only the case m = for which j(k) = 1 and time-reversal 
symmetry is possible with two component spinors. 
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Applying T to the plane-wave solutions on a bond (O) we find 

1 



-l fi 



e ifcx + i/t ' ' 
' 1 i 



\kx 



(4.21) 



Time-reversing spinors on all the bonds meeting at the vertex i defines new vectors of 
outgoing and incoming waves ~jx r , fj, T , 



'-Vli 







T(«i-u>j) 



-In 







I(t)j— Wi) 



For the vertex transition matrix to be time-reversal invariant we require 



(4.22) 



(4.23) 



Using the definition of the vertex transition matrix, ~[t = T^*jZ for unitary, equation 
( ggjj ) implies 



(T (i ))T 



'■Wi 





I 





(Vi-Wi) 







I 





(l>i-tOi) 



(4.24) 



Equation ( [4.24| ) is the condition that all vertex transition matrices must satisfy in order 
that the system possesses time-reversal symmetry. 

Splitting TW into four blocks equation ( 4.24 ) is equivalent to 



I'M 



Ti T 2 
T 3 T 4 



T 



-Ti 
T T 2 



TI 
-TI 



(4.25) 



The components T^ c of the transition matrix are thus either symmetric or antisym- 
metric according to the alignment of the bonds at the vertex, in the sense that: 



Tbc — T 



T 



be 



cb 



-cb 



b 

->- 



b 

-<r- 



or 



b 




This raises a number of problems. Back scattering is not possible in this scheme as is 
identically zero (returning to the same bond always requires antisymmetry). Consequently 
it is not possible to quantise graphs containing vertices with valency one. It is also clear 
that the components of TW cannot be invariant under a permutation of the bonds. But 
it is this natural physical property imposed by Kottos and Smilansky [Tj| [T7| that allows 
them to derive a simple form of the trace formula for the Schrodinger operator on a graph. 
In fact we may ask if any graphs exist for which the vertex transition matrices are both 
unitary and satisfy the time reversal condition ( (4.24j ). There are no examples of TW for 
vertices with valency one and only a trivial solution for one incoming and one outgoing 
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bond, Vi = 2. With Vi = 3 there are again no solutions. The simplest nontrivial example 
of a transition matrix is for a vertex with valency four. In this case, with two incoming 
and two outgoing bonds, one possible choice of a time-reversal invariant transition matrix 
is 



-1 
1 





1 

7! 



1 
o 
o 
i 



i 
o 
o 



(4.26) 



o\ 
1 
1 

-1 0/ 

Clearly this is both a unitary matrix and has the symmetry and antisymmetry properties 
prescribed by ( 4.25 ). Physically this would represent strange boundary conditions as a 
spinor reaching the vertex would be prevented not only from back scattering but also from 
scattering down another of the bonds. However it is still reasonable to ask how closely the 
spectrum of a graph quantised with such a vertex transition matrix will follow the GSE 
prediction. In the following subsection we investigate one example of a binary graph (i.e. 
where all vertices have two incoming and two outgoing bonds). 



4.3 Bond scattering matrix 

To find the energy levels we must first determine the bond scattering matrix § for the whole 
graph from the transition matrices at the vertices. Let us define //^ to be the coefficient 
of the plane wave on the bond (ij) traveling in the direction % — > j. Then 

is £(y) in 

our previous notation. The bond scattering matrix § = (S(ij)(i m )) is defined by 

S (mm )(k)=5 mt T$ ilm) e ikL ^ . (4.27) 

To scatter from the spinor with coefficient /i^ m ^ to that with fj,^ the bonds (Im) and 
(ij) must be connected at m. The transition amplitude is then defined by the transition 
matrix T^ m ^ at m. Before the transition the spinor collects a phase propagating along 



(Im). Equation ( f4.27| ) is equivalent to the description of the scattering matrix as a product 
E>(k) = TB(fc) where B>(k) is a diagonal matrix of phases and T is a matrix of transition 
amplitudes for the whole graph. As there are two coefficients for the spinors on each bond 
§ is a 2B x 2B matrix. The unitarity of the vertex transition matrices and the diagonal 
matrix of phases D(/c) ensures that §(/c) is also unitary. Note that the vertex transition 
matrices T^ m ^ and hence the bond S-matrix still depend on the directions assigned to the 
bonds. 

The bond S-matrix E>(k) acts on the vector of coefficients /i, which is composed of the 
coefficients ~jZ and ^7 at every vertex and therefore defines the plane wave on the graph. 
Energy eigenf unctions correspond to vectors \x where 

S(Jfe)/i = /x . (4.28) 

The energy eigenvalues then correspond to the values of k for which 

\hB -S(fc)| = . (4.29) 
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Figure 2: The integrated nearest neighbour spacing distribution (left) and a histogram of 
the spacing distribution (right) for a fully connected Dirac pentagon, with 35753 levels. 

Figure |2| shows the nearest neighbour level spacing statistics for a fully connected pen- 
tagon calculated using the vertex transition matrices ( 4.26Q . Directions were assigned to 



the bonds to produce two incoming and two outgoing bonds at each vertex. The assign- 
ment of TW to the vertices is still not unique as the elements of TW also vary between 
bonds with the same direction. The numerics confirm that, even for this small graph 
with unusual boundary conditions, the energy level statistics correspond well to those of 
random matrices from the GSE. We now turn to consider how the quantisation of more 
general graphs can be realised using two component spinors and whether there exists a 
trace formula for the Dirac operator. 



4.4 Paired bonds 

We expect it should be possible to put a Dirac operator on any (topological) graph while 
preserving time-reversal symmetry; and the assignment of elements in the transition ma- 
trix should follow a general scheme independent of the particular bonds to which they 
correspond. We saw that difficulties arise due to the fact that, as opposed to the case of a 
Schrodinger operator, the Dirac operator is a first order differential operator, thus breaking 
the invariance under x i— ► —x. In this section we show that a quantisation in terms of a 
Dirac operator can be achieved for two component spinors by replacing each bond of the 
classical (topological) graph with a pair of metric bonds one running in each direction. 

We maintain the realisation of a two component Dirac operator as described above, 
however with the constraint that each vertex i has valency 2vi with pairs of directed bonds 
meeting at i as indicated in figure [3| The bonds in a pair replace a single bond of the given 
topological graph and hence are assigned the same length but run in opposite directions. 
This restores the physical symmetry of the problem as all the original bonds are treated 
equivalently. Putting plane- wave solutions, equation ( |4.3| ), on the bonds we can write an 
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Figure 3: An isolated vertex with paired bonds. 



eigenfunction on a pair b of bonds as 



Akx b , 6 / 1 \ n ikx h 

e ~r Hp 



(4.30) 



-ikxi. 



where x Vi+ b has been replaced by L^ — Xb and the extra phases absorbed into the coefficients. 
The coefficients have been labeled a or /3 according to the type of spinor, characterised by 
a lower component i or — i, rather than the bond on which the spinor travels. Now let us 
rearrange the coefficients in the vectors 7^ and so that a coefficients are always listed 
first, 

( &\ 



( Hi\ 

Hp 



Hp 



(4.31) 



V : 7 



The self-adjoint boundary conditions still define a unitary vertex transition matrix, 

~Ji = T (i) V and (T^)- 1 = (T^) f . (4.32) 

Applying the time-reversal operator Q4.20 ) to the wavefunctions ( |4.30|) we find time- 
reversed vectors 7* r> ^r> 



/ Hp \ 



-Hi 



\ 



( h\_\ 
-Hi 



J 



(4.33) 



Time-reversal invariance again implies 



Mr 



pj(i) 



(4.34) 
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Using the definition of the transition matrix (|4.32|) we obtain the following condition for 
time-reversal invariance, 



(4.35) 





( J" 1 






( J 


\ 
















V 











where 



J :-- 



1 
-1 



(4.36) 



The form of equation ([4.35 ) is similar to, but not the same as, the definition of a symplectic 
matrix. 

The vertex transition matrix can be divided into 2x2 blocks (T^~)) bc which relate 
incoming spinors on the pair of bonds b to outgoing spinors on the bond pair c. Let 



^rjp(i)\bc 

then time-reversal symmetry implies 
or equivalently, 



r s 
V w 



w —s 
—v r 



(T^) cb = |(T^) bc | ((T (i )) 6c ) 1 



(4.37) 

(4.38) 
(4.39) 



This suggests both a method of constructing time-reversal symmetric transition matrices 
and the form of permutation invariance among the bond pairs to be required. Matrices of 
this form can be written 



( Ui 

V 



\ 



Vi ) 



X 



1 
1 



\ 



u 



\ 



Vi j 



(4.40) 



where X = (xb c ) is a symmetric unitary matrix of dimension Vi and Uj G SU(2). \xb c \ is the 
probability that any incoming state on the classical bond c scatters to any outgoing state on 
the classical bond b. We may now choose to implement a physical permutation symmetry 
between the bonds of the classical graph by requiring the matrix X to be invariant under 
permutations a, 

X = aXa- 1 . (4.41) 
Together with the condition that X is unitary this implies 



X = e w 



( 9-1 



V 9 



Q 



\ 



with q :- 



1 + e* 



(4.42) 
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Unsurprisingly this is similar to the vertex transition matrix for a Schrodinger operator 
fl6| , |17| which is also symmetric unitary and permutation invariant. Here the coefficient 
lu, which we can think of parameterising the boundary conditions, is independent of k. 
The extra phase e l6> is not interesting as it corresponds only to an extra k independent 
phase factor. Adjusting the parameter ui we see that the Dirac equation also allows both 
Neumann like (u = 0) and Dirichlet like (u = tc) boundary conditions. An element u^u~ x 
of SU(2) defines a 'spinor rotation' from the pair of linearly independent two component 
spinors incoming on bond c to the outgoing pair on bond b. 

Equations ( 4.40 ) and (|4.42|) define general time-reversal invariant vertex transition ma- 
trices for two component spinors. The transition amplitudes have a form of permutation 
invariance ( (4.41|) and in section |6.1| we see that this allows a simple form for the trace for- 
mula to be derived. With the vertex transition matrix the bond scattering matrix for the 
whole graph can be defined as in equation ( 4.27Q with unpaired bonds. The quantisation 
condition ( |4.29| ) is also unchanged, apart from the fact that due to doubling of the number 
of bonds E>(k) now is a 4B x 4B matrix. 



5 The Dirac operator for four component spinors 

Before looking at examples for the two component spinors on paired bonds we will compare 
with the alternative construction of the Dirac operator acting on four component spinors. 
Therefore, we now consider again directed graphs with B unpaired bonds. In this case the 
time independent Dirac equation reads 



•fc d ^ , 2 a I 

—inca— V mc pip 

ax 



(5.1) 



where we choose a and (3 as in 

/ 



a 



\ 





— i 




h 

-i\ 







p 



( 1 









-1 











(5.2) 



1 J 



We could alternatively have chosen either of the other a matrices of the standard represen- 
tation for the Dirac operator in three dimensions or used another representation entirely. 
The choice here simplifies the working. 

Eigenspinors of T> are again plane waves. For positive energy they are of the form 

/ 1 \ / \ 

° e lkXb + n b J 



fe) 



o 

V i7(*0 J 

( 1 





\ 



ij(k) 


/ 



Akxh 



-ikxi, 



J 

o \ 



(5.3) 



1 

vy(k) 




-ikxi, 
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with k > 0. The dispersion relation and j(k) are as in ( |4.4j ). These four plane-wave 
solutions differ from the two component case although they are labeled with coefficients 
li h a etc. in analogy with the two component solutions on pairs of bonds. The wavefunction 
ip b (xb) is a function on a single directed bond of a metric graph. We notice that, in contrast 
to the pairs of two component wavefunctions ( |4.3C|) , ip (#&) is not invariant under a change 
in the bond direction x b \— > —x b . 

5.1 Self- adjoint realisations on graphs 

The construction of self-adjoint realisations of the Dirac operator on the graph proceeds 
in complete analogy to the case of two component spinors. We therefore must look for 
maximal closed subspaces of 

B 

0W 2> i([O,L 6 ])®C 4 (5.4) 

6=1 

on which the skew-hermitian form 

n{M) = w,v>-<0,zty> (5.5) 

vanishes. In order to convert the problem to the boundary values of the spinors we integrate 
by parts and express Q as a complex symplectic form, 

on C 8B . Here we have defined 

V> + := tyi(0), • • • , (0), ^(0), . . . , Vf(0), 

^\{L X ), . . . , ^ ^(Lx), . . . , ^{L B )f , 

V>~ := (-Vl(O), . . • , -<(0), ^(0), • • • , Vff(0), 

^(Li), . . . , tpf(Ls), -^{Lrl . . . , -<l>f(L B )) T • 

The map from the spinors on the graph to the vectors V + an d ip~ mixes the first and 
second components of the spinors in and the third and fourth ones in . Identifying 
maximally isotropic subspaces of C SB with respect to Q, equation (|5.6| ), the argument 
proceeds as for the two component spinors. Such a linear subspace is defined by 

A^ + +B^"=0, (5.8) 

with 4B x AB matrices A and B, and is maximal isotropic if rank (A, B) = 4B and AB^ = 
BA'. Under such conditions the Dirac operator on (|5.4j ) with boundary conditions as 
prescribed by (|5.8|) is self-adjoint. 
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For eigenf unctions of the Dirac operator, i.e. plane-wave solutions ( |5.3| ) on each of 
the bonds, we define vectors of incoming and outgoing coefficients at a single vertex, for 
example figure [I], 



/ 



Ma \ 
/4 



LL Wi 
Ha 

Wi 

Ha c 



( 



r.wi 

Ha 



\ 



H'a e 

e ife ^i+i 



e ikLv i 



(5.9) 



/ 



Then using boundary conditions defined at the vertex by matrices Aw and Bw we find as 
in the two component case, 



Jt = -(A w - i 7 (A;)B (i) )- 1 (A w + i 7 (A;)B (i) ) /I . 
The vertex transition matrix 

T« : = -(A w - i 7 (A;)B w )- 1 (A (i) + i 7 (fc)B w ) 



(5.10) 



(5.11) 



is unitary due to the condition A^BW — ]gW^Wt that guarantees current conservation at 
the vertex. 



5.2 Time-reversal symmetry 



We now turn to see what time-reversal invariance implies for the vertex transition matrix. 
The time-reversal operator in the standard representation is 



r 



( o 


— i 








\ 


i 






















— i 









i 








(5.12) 



This commutes with the Hamiltonian on the bonds so there is no requirement that the 
mass be zero with four component spinors. Applying T to the wavefunction at the vertex 
i we define new vectors from the coefficients of outgoing and incoming waves, ~jtq-, "/^V, 



/ J 



\ 



J ) 



( J 



Mr 



\ 



J ) 



(5.13) 
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For time-reversal invariance we require 

7? r = T (i) ^ r . (5.14) 

Using the unitarity of time-reversal symmetry ( 5.14j ) is equivalent to demanding the 
transition matrix satisfy 



(5.15) 



This is the same as condition ( |4.35 ) derived for two component spinors on paired bonds. 
The general form of transition matrices ( |4.40| ) proposed for the spinors on paired bonds 
will hence also hold for four component spinors on a directed graph. Furthermore with 
four component spinors the mass is no longer required to be zero. Instead any which 
can be constructed from the boundary conditions according to equation ( |5.11| ) and which 
is also time-reversal invariant is available. 
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6 Energy level statistics 



To verify the conjecture of Bohigas, Giannoni and Schmit in the case of the time-reversal 
invariant Dirac operator on a graph we examine the energy level statistics through both 
numerical calculations and directly via the trace formula. 

From now on we consider only boundary conditions with the form of general vertex 
transition matrices defined in equations (f4.40|) and ( |4.42|) . This still leaves a range of 
boundary conditions parameterised by uo. Of particular interest are the Neumann like 
boundary conditions, uj = 0. These are the boundary conditions most often studied for 
the Schrodinger operator on the graph. The vertex transition matrix for Neumann 
boundary conditions is generated by the matrices 

/ 1 -1 ... \ 
-1 ... 















1 









( 



V 



J 



where is a block diagonal matrix of elements of SU(2), 



V 



U Vi J 



(6.1) 



(6.2) 



(6.3) 
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Substituting these boundary conditions into equation ( |5.11| ) and using some algebra we 
obtain the vertex transition matrix 



/ 3- 



V 



- \ 



/ 




(6.4) 



TW is independent of the function 7(fc) determined by the mass. With this form of the 
Neumann boundary conditions there is no need to specify the mass or whether the system 
is represented by paired two component spinors or four component spinors on a directed 
graph. 
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Figure 4: The integrated nearest neighbour spacing distribution and a histogram of the 
spacing distribution for a Dirac operator with Neumann boundary conditions on a fully 
connected square, 24234 levels. 

Figure |] shows the nearest neighbour spacing statistics for a fully connected square with 
Neumann boundary conditions, the vertex transition matrix defined in equation (|6.4j) . The 
three elements w& G SU(2) at a vertex were chosen as 

u b = exp(i0 b o- n J , (6.5) 

where a nb is a Pauli matrix. The parameters 9b were then selected randomly at each vertex. 
To obtain good agreement with the symplectic random matrix ensemble the parameters 
were also picked such that the elements Uf, generated at the vertices were all sufficiently 
different. Deviations from the GSE behaviour become apparent if the elements of SU(2) 
at one vertex are similar to those at another, i.e. 9\ ~ 9™. 

Introducing a magnetic vector potential Af, on the bonds breaks time-reversal symmetry. 
The transition elements for the ^-matrix remain the same but a plane wave propagating 
down a bond (i — » j) picks up an extra phase exp(L4(y)L(^)) where we know that Lu^ = 
L(jj) but A(ji) = —A(ij). Figure || shows that for the Dirac operator on a fully connected 
square breaking the time-reversal invariance in this way produces energy level statistics 
like those of random GUE matrices as expected. For the Dirac operator we used Neumann 
boundary conditions and random elements of SU(2) chosen as previously. 
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Figure 5: The integrated nearest neighbour spacing distribution and a histogram of the 
spacing distribution for a Dirac operator with Neumann boundary conditions on a fully 
connected square where time reversal symmetry has been broken with a magnetic vector 
potential, 25503 levels. 

6.1 Trace formula 

A trace formula for the Laplacian on a graph was first produced by Roth ||18|| . We derive 
the trace formula for the density of states of a system with Neumann boundary conditions, 
following the approach developed by Kottos and Smilansky Jl6|, [17| for the Schrodinger 



operator. The vertex transition matrices for such a system were defined in equation 
The main tool in setting up the trace formula is the bond scattering matrix S(fc) defined 
in equation ( |4.27| ). 

Consider the density of states for the wave-number k. It is an infinite series of delta 
functions located at k = k n corresponding to eigenvalues E n of the Dirac operator. The 
eigenvalues are counted according to their multiplicity. Wave-numbers k n are values of k 
for which the function 

C(A:):=|l4B-S(A;)| (6.6) 
is zero. Diagonalising §(fc) we can write ((k) in terms of the eigenvalues e 1 ^^ of S(k), 

\hs - S(*)| = fl(l - eW*>) = 2 4B |§(^ ffsin (^f) . (6.7) 

C(k) is complex with the phase contained in the term |S(fc)|i. Multiplying by |S(fc)|~^ we 
define a real function whose zeros are also the zeros of C(^)- Then the density of states 
can be written 

d(k) = y^5(k-k n ) = limlm— log(|l4B-S(Jfe + ie)||S(A: + ie)|-3) . (6.8) 

Z — ' 71 e^O dk 

n=l 

The expression (p.8[) factorises into a sum of two terms. The smooth part of the spectral 
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density is 



d sm th(k) := --limIm-^-log(|§(£; + ie)| *) 
2L 

5 

7T 



(6.9) 



where L is the total length of the graph, 

B 
6=1 

(Note that in our definition of L the length of each bond of the classical graph is counted 
only once although the S-matrix can be considered either as scattering between two com- 
ponent spinors on paired bonds or four component spinors on a directed graph.) 
The oscillating part of the spectral density is 

d osc (k) :=--UmIm^-log(|l4B-S(A; + ie)|) . (6.10) 
The logarithm can be expanded as a sum of powers of the trace of S(k), 

CO 

log(|I 4B - S(fc)|) = - V - trS n (A;) . (6.11) 

L — ' n 

n=l 

As in the case of the Schrodinger operator on a graph we can now write the trace of powers 
of S(k) as a sum over periodic orbits on the graph. However, for the Dirac operator, the 
transition amplitudes between classical bonds are specified by the 2x2 matrices T 6c , 

tr§ n (A;) = tr ^ e ^L blTbl b 2e ikL b2 jb 2 b 3 _ _ _ ^kL^jb^ _ ^ ^ 

bi...b„ 

Thus tr E> n (k) ^ implies that the sequence (bib 2 . . . b n ) of bonds labels a periodic orbit on 
the graph. For a single periodic orbit p which passes through n vertices, 

tr(T M2 . . . T Kbl ) = A p tr(d p ) e i7r ^ . (6.13) 

In this formula A p is a (positive) stability factor, the product of the absolute values of the 
transition amplitudes for spinor pairs round the orbit, 

v=§H)ng), <iu, 
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where fj, p is the number of cases of back scattering from vertices with v s > 2 and n = fi p +v p . 
An overall sign is collected in the factor e 17TtJ-p and d p is the product of the elements of SU(2) 
that rotate between the two types of spinors on the periodic orbit, 

d p :=u blb2 u b2b3 ...u Kbl . (6.15) 

The spin rotations u bibi+1 := Ub i+1 u^. are composed of SU(2)-blocks comprising the matrix 
, see equation ( |6.3D , attached to the ith vertex visited along the periodic orbit. 
Combining these results, 

trS n (A;) = A p e ^ p tr(dp) , (6.16) 

peP„ 

where the sum is over the set P n of periodic orbits of n bonds and l p is the metric length 
of the orbit p. Substituting equations ( |6.16|) and ( |6.11| ) into equation ( |6.10| ) we find the 



oscillating part of the density of states. Adding the result for the smooth part ( |6.9| ) we 
obtain 

or i / 

d(k) = — + - V -^A v e m ^ tr(<L) cos(ifeL) . (6.17) 

7T 7T r p 



The sum is over all periodic orbits. If an orbit p consists of r p repetitions of a shorter 
orbit then only the length of the primitive orbit l p /r p is included in the formula. This 
trace formula provides an exact relation for the density of states of the Dirac operator on 
a graph with Neumann boundary conditions as a sum over the classical periodic orbits. 
Equation (|6.17| ) for the density of states can be compared both to the form derived for 



the Schrodinger operator on a graph |jl8| , |I6| , |l7f and to previous semiclassical results for 
systems with spin ]7|, [8], f[4j . Most terms in our density of states are the same as those 
derived by Kottos and Smilansky for the Schrodinger operator [IB, 17 1. The additional 
term for the Dirac operator is the trace of d p . This can be thought of as an additional 
weight factor which depends on the transformation of the spinor pairs around the orbit. 

Comparing our trace formula to the general semiclassical form for the Dirac operator 
]7J, ||, we see that here the Dirac operator also produces a trace over an element of SU(2) 
associated with the change in spin around the orbit. This weight factor can also be inter- 
preted in terms of an effective rotation angle of a classical spin vector transported along 
the periodic orbit according to the equation of Thomas precession. That the results agree 
is not unexpected. The Dirac operator on a graph however does not require a semiclassical 
approximation to derive the density of states. Finally we also note the similarity with the 
trace formula for a different system, the quantum cat map with spin, derived by Keppeler, 



Marklof and Mezzadri, [TJj]. In this case the quantum map is an equivalent to the Pauli 
operator with spin 1/2. In their semiclassical trace formula the presence of spin also ap- 
pears as a trace of an SU(2) matrix which defines the spin transport around a periodic 
orbit. 
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6.2 Form factor 



From the trace formula we derive the spectral two-point form factor for the Dirac operator 
on quantum graphs. The form factor is an energy level statistic that is widely studied 
for classically chaotic quantum systems and in particular for the Schrodinger operator on 
graphs, see fIB], [17|, |, |], ||. We see that our form factor agrees with the general semiclassical 



result for systems with spin ||, and making the diagonal approximation the form factor 
for low r is consistent with the GSE prediction from random matrix theory. 

We derive the form factor for the energy level spacings. This is the Fourier transform 
of the two-point correlation function itself defined as 

R 2 (x):=( — ) lim - / d(k) d(k —) dk — 1 . (6.18) 

\2LJ a^oo A Jq V L / 

In defining the two-point correlation function we have rescaled the spectrum dividing by the 
mean level spacing 2L/n, see equation ( |6.9| ). For systems with half- integer spin and time- 
reversal invariance the eigenvalues k n come in pairs (Kramers' degeneracy). We consider 
the spectrum x n where this degeneracy has been lifted, each pair of eigenvalues being 
replaced with a single representative. The mean spacing of the x n distribution is therefore 
L/n. Substituting the trace formula for the density of states (|6.17 ) into the definition of 
R2(x) and carrying out the integral we find 

The two-point form factor K 2 {r) is the Fourier transform of the two-point correlation 
function, 



K 2 (t):= R 2 (x)e 2niXT dx , (6.20) 

J — oo 

which is even in r. Taking the Fourier transform of R 2 (x) term-wise we obtain 



K 2 (t) 



1^2 E AA e^ + ^ tr(d p ) tr(d q ) 5 (t - ^) 5 lp , k (6.21) 



4(2L) M . plq 



for r positive. 

The eigenvalue distribution of random matrices averaged over the GSE has a two-point 
form factor 

^ (T)= MH-J|r|log|l-|r|| forjrj^ 
z y ' \ 1 for \t\ > 2 v ' 

Making similar assumptions to those often used in the Schrodinger case we will see that 
K2{t) ~ M/2 for |r| — ► 0. This agrees with the random matrix theory prediction of 
equation ( |6.22| ). 

The analysis of the form factor based on the trace formula requires us to take a semi- 
classical limit. The latter can be characterised by going to some regime with increasing 
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mean spectral density. According to the relation ( |6.9| ) this corresponds to increasing the 
total length L of the bonds. We realise this limit in passing to graphs with an increasing 
number B of bonds while keeping the mean bond length L := L/B constant. As r — > 
the form factor should be dominated by the diagonal form factor introduced by Berry ||, 

^^iS 5 ^^) {A>(d p )) 2 + A p A F tr(d p )tr(rf F )}^r-^ , (6.23) 

since in this limit only correlations for an orbit p with itself and with its time-reversed 
partner p are relevant. Combined with the semi classical limit we now have to consider 
L — > oo, t — > with Lt — > oo. 

Comparing p and p we see that A p = A p and dp = d~ x . For elements of SU(2) tr(rf _1 ) = 
tr(d), 

^W = ^Eg) 2 4W4))^(-^) • (6.24) 

In the limit Lt — > oo it is the long orbits which dominate the sum and for these the 
proportion which are repetitions of shorter orbits tends to zero so we ignore periodic orbits 
with r p 7^ 1. Breaking the sum on periodic orbits into a sum on sets P n of orbits of n 
bonds and using l p m nL, 

Kf a9 ir) « \ J2 (^) 2 5 (r - A) ^ ^ (tr(^)) 2 . (6.25) 

We can regard both A p and d p as functions of p. The elements of SU(2) in d p were 
selected randomly, according to Haar measure, so A p and d p are uncorrelated. If the 
number of periodic orbits of length n is \P n \ then for large n 

1 n| pePn V n| pe-Pn / V n| pePn / 

We now consider the two averages separately. The sum over P n of A 2 p can be evaluated in 
the limit n — > oo using the classical ergodicity of the graph. Classical motion on the graph 
is determined by the doubly stochastic matrix M = (M bc ), with the magnitude squared of 
the transition amplitudes between the bonds as entries, 

M bc := \X bc \ 2 . (6.27) 

M defines a Markov chain. For a connected graph a given pair (b, c) of bonds can be linked 
by a path of q bonds. Thus {M q ) bc > such that M is an irreducible matrix and the 
Markov chain is ergodic. If, in addition, no eigenvalue other than one lies on the unit circle 
it is mixing which implies 

lim (M") 6c = -L . (6.28) 



2B 
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The trace of M n therefore approaches unity, but the trace can also be expressed as a sum 
over the periodic orbits, 

tiM n = n^A 2 p . (6.29) 

The factor n counts the cyclic permutations of the orbit. Thus in the limit of long orbits 
we may replace ^2 pePn A 2 with 1/n if the Markov process is mixing. It should be noted 
that while we have used the mixing property of the classical motion on the graph without 
spin (defined by the matrices X) it is equivalent to the same property of a 4B dimensional 
matrix M defined from T. In this case Mf, c = \Tt> c \ 2 and each spinor is treated separately. 
We defined the mixing property of the classical graph using X to make the connection 
with the amplitudes A p clear. 

In the limit of n — > oo we replace the average over P n of (tr(<i p )) 2 with the integral over 
SU(2) with Haar measure. Evaluating the integral directly we obtain 



lim 

n— >oo 



Pn\ 



pEPn 



(tr(u)) 2 du 



(6.30) 



SU(2) 



Substituting these results into the diagonal approximation for the form factor fl6.25|) we 
find 



K diag { 



T ~ 



-H 1 

9 L^i Of 



n 



IB IB 



5 



n \ 
T -2B) 



1 



T 



(6.31) 



in the combined limit L — > oo, r — > with Lt — > oo. This is in agreement with the GSE 
prediction for small \t\. The validity of the assumptions made in the argument is of less 
consequence than that they are of the same type used in other semiclassical approximations, 
particularly for the Schrodinger operator on a graph. The key additional assumption made 
to include the spin evolution was equation ( |6.26| ), taking the elements d p of SU(2) to be 
uncorrelated with the weights A p of the orbit. 
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Figure 6: Plots of an averaged spectral form factor from the fully connected square and 
the GSE prediction. The inset shows the result for the form factor without averaging. 
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Figure |^ compares a numerical calculation of the form factor with the GSE prediction 
( |6.22j ). As the form factor takes the form of a distribution, see equation ( |6.21| ), meaningful 
results are only obtained when it is averaged with a test function. For simplicity we use 
a function which is one in the region r ± 0.1 and zero outside. The theoretical prediction 
for the GSE form factor can also be averaged in the same way. As the computation of the 
form factor leads to a delta function at r = in the limit of the number of levels going to 
infinity the first few values of r are omitted in the numerical calculation. 



7 Gedanken experiment 



It is common when studying graphs, in this context, to describe them as an idealised model 
of a "typical" quantum system. When compared to quantum maps or systems of constant 
negative curvature the case is certainly strong. However in introducing the Dirac operator 
to graphs we may have inadvertently weakened the argument and that is what we seek to 
correct here. 




Figure 7: A graph realised as a network of wires embedded in 3d. 

Let us think of the classical graph that we wish to quantise as a network or web of 
wires set up in our ideal laboratory, see figure [7[ To quantise this system using the Dirac 
operator we must introduce pairs of spinors propagating in both directions along the wires. 
At the vertices changes in the composition of the spinor pairs are described by elements 
Ubc € SU(2), c labels the incoming classical bond and b the outgoing bond. These elements 
Ub c necessarily depend on the pair of wires under consideration, w fec ^ w rfe . This is not a 
feature of the usual Schrodinger quantisation [16, 17] where all wires meeting at a vertex 



are treated equivalently. It is therefore fair to ask whether our model web could really 
distinguish bonds in such a way. 

One simple geometric argument connecting the matrices to our experimental net- 
work is provided by the double covering map / : SU(2) — ► SO (3). Let R(9b c ) G SO (3) 
describe the rotation from bond c to b at a vertex. Setting u 6c = /~ 1 (i?(6 l bc )), with some 
choice of the branch / _1 (SO(3)), we define a set of elements of SU(2) from the architecture 
of the network. This construction is time-reversal symmetric as u c b = u^. (To write this 
in the form used for the vertex boundary conditions (|6.4j ) choose a reference direction at 
each vertex and take the v elements u c to be f~ 1 (R(6 c )) where R(6 C ) rotates from the 
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reference direction to bond c.) We have not attempted to provide a physical justification 
for this transformation of spinor pairs at the vertices. Our formulation rather shows that 
the extent to which bonds must be distinguished when quantising a graph with the Dirac 
operator remains consistent with the picture of the graph as a simple ideal quantum sys- 
tem. There is sufficient geometrical information in the structure of a graph embedded in 
three dimensions to determine rotations of the spinors at the vertices. 
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